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Abstract Consider the control system (E) given by x — x(f + ug), where x £ SO(3), \u\ < 1 and 
/,j 6 so(3) define two perpendicular left-invariant vector fields normalized so that |j/|j = cos(a) 
and \\g\\ = sin(a), a e]0,7r/4[. In this paper, we provide an upper bound and a lower bound for 
N(a), the maximum number of switchings for time-optimal trajectories of (E). More precisely, we 
show that Ns(ct) < N(a) < Ns(ct) + 4, where Ns(a) is a suitable integer function of a such that 

N s (a) a * o 7r/(4a). The result is obtained by studying the time optimal synthesis of a projected 
control problem on MP 2 , where the projection is defined by an appropriate Hopf fibration. Finally 
we study the projected control problem on the unit sphere S 2 . It exhibits interesting features which 
will be partly rigorously derived and partially described by numerical simulations. 

1 Introduction 

Let (E) be the control system given by: 

x = x(f + ug), (1) 

where x G 50(3), |w| < 1 and f,g € so(3) give rise to two non zero perpendicular left-invariant vector fields 
on SO(3). In this paper, we consider the following problem: given any pair of points 21,2:2 of SO(3), find a 
trajectory of (1) steering x\ to x-i in minimum time. That issue is known as the problem of determining the 
time optimal synthesis (TOS) for (E). The strategy to determine a TOS usually consists in two steps: 

1. Reduction procedure: it is based on the Pontryagin Maximum, Principle (PMP) which is a first-order 

necessary condition for optimality. Roughly speaking, the PMP reduces the candidates for time optimality 
to the so called extremals, which are solutions of a pseudo-Hamiltonian system. This reduction procedure 
may be refined using higher order conditions, such as Clebsch-Legendre conditions, higher-order maximum 
principle, envelopes, conjugate points, index theory (cf. for instance [2, 3, 5, 12, 18, 19, 21, 26, 27, 28, 
30, 32, 33, 34, 35]); 

2. Selection procedure: it consists of selecting the time optimal trajectories among the extremals that passed 

the test of Step 1. (see for instance [8, 12, 16, 25]). 

Step 1. is already non trivial and, in general, the second one is extremely difficult: if the state space is 
two-dimensional, the problem of determining the TOS for single-input control systems is now well-understood 
[9, 11, 12, 14, 15, 23, 24, 32, 33]. However, for higher dimensions, very few examples of complete TOS for a non 
linear control system are available (see for instance [29]). Intermediate issues were thus deeply investigated: 
determining estimates for the number of switchings of optimal trajectories, describing the local structure of 
optimal trajectories, finding families of trajectories sufficient for optimality, cf. [6, 13, 19, 22, 26, 28, 36], etc. 

For the control system (E), we normalize the two perpendicular vector fields induced by / and g in such a 
way that ||/|| = cos(a), ||g|| = sin(a), with a s]0,7r/2[ (for the precise meaning of "perpendicular" and of the 
symbol ||.||, we refer to Section 2.1). Defining X + := f + g and X_ := / — g, we have \\X + \\ = \\X-\\ = 1 and 
a is the angle between / and X + . 

By a standard argument (see Section 2 below), one can show that every time optimal trajectory is a finite 
concatenation of bang arcs (i.e. u = ±1) or singular arcs (u = 0) and thus, the Fuller phenomenon (i.e. 
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existence of a trajectory of a control system joining two points in (finite) minimum time, with an infinite 
number of switchings, cf. [20, 37]) never occurs. (A switching time - or simply a switching - along an extremal 
is a time to so that the control u is not constant in any open neighborhood of to-) Moreover, one can easily 
show that, the supremum N(a) of the number of switchings over all time optimal trajectories of (£), is finite. 
By using the index theory developed by Agrachev, it is proved in [3] that 

N(a) < N A := f-1 , (2) 
La J 

where [•] stands for the integer part. That result was not only an indirect indication that N(a) would tend to 
oo as a tends to zero, but it also provided a hint on the asymptotic of N(a) as a tends to zero. 

A related line of work regards the study of the distributional version of (£), which is the driftless control 
system given by x = x{uifi+U2f2), |ui|,|u2| < land/1,/2 € so(3) linearly independent. Indeed, assuming that 
1 1 /ill = 1 1/2 1|, Sussmann and Tang ([36]) showed that time optimal trajectories have at most four switchings and 
they provided a finitely parametrized family of trajectories sufficient for optimality. That result was extended 
to the general case (/1 and f 2 just linearly independent, cf [13]): time optimal trajectories have at most five 
switchings. For both works, the elimination from optimality of extremals with respectively five or six bangs 
relics on the envelope theory developed in the context of control theory by Sussmann (cf. [35]). 

At the light of the previous results, there was strong evidence for two radically situations as a tends to zero: 
for (£), N(a) is expected to go to infinity, as for the distributional control system, there exists a universal 
bound on the number of switchings. The main result of the present paper confirms that difference, i.e. N(a) 
tends to 00 as a tends to zero. More precisely, we complete the inequality (2) as follows: 

Theorem 1 Let (£) be the control system defined in (1) with f,g perpendicular so that ||/|| = cos(a) and 
\\g\\ = sin(a), a e]0,7r/4[. Then, if N(a) is the maximum number of switchings along a time-optimal trajectory 
of (£), we have 



N s (a) < N(a) < N s (a) +4, where N s (a) := 2 — - 2 — . (3) 

-8a J L L8aJ 4a. 

The above theorem improves (2) in two ways: i) for a small, it (essentially) divides the upper bound of N(a) 
by four with respect to (2); ii) it provides a lower bound of N(a) differing from the upper bound by a constant. 

The lower bound is in fact our main contribution and, to get it, one must prove the existence of time optimal 
trajectories of (X) admitting at least a number of switchings equal to that lower bound. Our strategy consists 
of projecting the control problem onto another (E)s defined next. First, let MP 2 be the two-dimensional real 
projective space (i.e. the two-dimensional manifold made of the directions of M 3 ) and fix a point x e SO(3). 
Consider the Hopf fibration II : SO(3) — ► RP 2 defined by KerdH(xo) = Span{xof}, which means, roughly 
speaking, that II annihilates the drift term / at xq. Then, we project (S) by II and obtain a single- input 
S'0(3)-equivariant control system (S)g on RP 2 given by y = y{fs + ugs), with fs — dH(f) and gs — dU(g), 
that is locally controllable. We then consider the minimum time problem of connecting Il(a;o) to any other 
point of RP 2 . 

In fact, we study a slightly different time optimal problem by lifting (£)s to the unit sphere S 2 . By an 
abuse of notation, we still denote by (£)s the control system obtained in that way. Hence H(xo) is identified 
with the north pole and MP 2 is identified with NT_, the subset of the sphere made of the union of NI, the 
(open) top hemisphere of S 2 , together with half of the equator. 

The time optimal problem consists now of connecting, in minimum time, the north pole with any point of NI . 
Thanks to the suitable choice of the Hopf projection and since a belongs to the interval ]0, 7r/4[, all extremals 
of the projected problem are bang-bang (i.e. they are a finite concatenation of trajectories corresponding to 
controls +1 or —1). Let Ns(a) be the supremum of the number of switchings for time optimal trajectories of 
(X)s starting at the north pole and ending in NT (such trajectories of (£)s are actually entirely contained in 
NI , see Lemma 7). 

The use of the Hopf fibration n is motivated by two facts: first, every time optimal trajectory for the 
time optimal problem on (£)s staying in NT is the projection by n of a time optimal trajectory for the time 
optimal problem on (£) and thus, Ns{ot) < N(a). Taking full advantage of the theory developed in [12], we 
will actually compute exactly N$(a) as given in (3). Second, using the fact that the fiber above H(xq) is the 
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support of a singular arc (for this problem singular arcs are integral curves of the drift xf) , we show that every 
regular bang-bang trajectory with at least Ns(a) + 5 cannot be optimal and thus, the upper bound. 

It is then clear, by now, that the most delicate part of the argument relies on the exact determination of 
Ns(a). This is done by studying the time optimal synthesis, (TOS for short) for the time optimal problem on 
(S) s . Such a TOS is usually constructed, following the theory developed in [9, 11, 12, 14, 15, 23, 24, 32, 33], 
recursively on the number of extremals arcs, and by checking at each step whether they are optimal or not. For 
the problem on RP 2 , we are not able to complete all the steps of the above construction, which would imply 
as a byproduct the existence of the TOS. In particular, we cannot show the optimality of all the extremals 
(i.e. the trajectories candidate for time optimality), but, from their study, we can demonstrate enough partial 
results in order compute JVg(a) precisely and thus to conclude the proof of Theorem 1. 

The complete time optimal synthesis is then studied numerically (actually on the whole S 2 ) and is showed 
in the top of Fig. 11. In particular, due to the compactness of S 2 , one of the main issues is to understand 
the singularities developed by the minimum time wave front as it approaches to the south pole. We provide 
numerical simulations that describe the evolution of the extremal front. As a — > 0, these numerical simulations 
suggest the emergence of three cyclically alternating patterns of optimal synthesis, each of them depending on 
an arithmetic property of a. 

The balance of the paper is organized as follows: Section 2 collects basic facts relative to the time-optimal 
trajectories of (£); in Section 3, the Hopf fibration is described and the proof of Theorem 1 is provided, assuming 
some facts about the time-optimal synthesis of (E)si whose arguments are deferred to the next section. In 
particular we use the expression of Ns(a) and the relation between the length of interior bang arcs for the 
problem on M.P 2 . The construction of the time-optimal synthesis of (£)s is investigated in Section 4, where an 
exact computation of Ns(a) is established. We conclude the section with two remarks, the first one explaining 
the relation between the TOS on the sphere and the TOS of a controlled linear pendulum, the second one 
establishing a link with an optimal control problem on SU(2). Finally, in Section 5, we provide the results 
of the numerical simulations completing the study of the time optimal synthesis (in particular of the possible 
behaviors in a neighborhood of the south pole) and we propose some open problems stated as conjectures. 



2 Statement of the Problem and Properties of Optimal Trajectories 
2.1 Basic Facts 

In this paper, we consider the control (£) given by (1), where x E SO(3), \u\ < 1 and f,g E so(3). An 
admissible control u is a measurable function u : [a,b] — ► [—1,1], where a, b depend (in general) on u, cf. [18]. 
A trajectory 7 of (£) is an absolutely continuous curve 7 : J — > SO(3), where J = [a, b] is a compact segment 
of K such that there exists an admissible control u for which 7(f) = 7(i)(/ + u(t)g) holds a.e. in J. We then 
say that (j,u), defined as before, is an admissible pair for (£). 

Definition 1 A trajectory 7 of (£), defined on [a, b], is time optimal if, for every trajectory 7' of (£) defined 
on [a' , b'] with 7(0) = 7' (a') and 7(6) = ^'(b 1 ), we have b — a < b' — a' . 

The Lie algebra (so(3), [., .]) is isomorphic to the Lie algebra (M 3 , x), where x denotes the vector product. This 
isomorphism is realized by the map: 

<Pl ■ so(3) -> R 3 (4) 




and provides an inner product on so(3) given by < Zi,Z2 >'■=< 4>l(z\), 4>l{z2) >, where Zi,z% € so(3). The 
symbol < ., . >, in the right-hand side of the above equation, stands for the Euclidean inner product of R 3 . 
With this definition, it follows that: < Zi,Z2 >■= —\Tr{z\Z2). In other words, this scalar product is the 
opposite of the Killing form on so(3). In the following, ||z|| := ^/< z, z > and Id is the 3x3 identity matrix. 
We will sometimes consider the 2x2 matrix corresponding to the planar rotation of angle (3 and we use Rp to 
denote it. 
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In this paper, we will assume that / and g are perpendicular and normalized so that ||/|| = cos(a) and 
||g|| = sin(a), a G]0,7r/2[. Here, we adopt the following notation used throughout the paper, c a :— cos(a), 
c 2 := cos 2 (a), s a :— sin(a) and s 2 := sin 2 (a). We define h := [f,g] = fg — gf and 

X + :=f + g, X_:=f-g. 

Note that \\X e \\ = 1, with £ = +,—. For a vector field z G so(3), we use e tz to denote the flow of z, acting on 
the right, so that t 1— > pe tz G 50(3) is the integral curve of z starting at p at time 0. Since z is linear, we have 
e tz — X)^Lo riT- We use ac ^ to denote the operator w ^ [z,w], acting on vector fields. If z,w are vector 
fields, then e* &<dz (w) :— e tz we~ tz . The Lie bracket relations between /, g, h are 

[f,g] = h, [g,h} = s 2 J, lh,f} = c 2 a g. 

From them, one deduces the following classical relations that will be useful later: 



e tadx '(/) 


= (c 2 a + 3 2 a cos(t))/ + ec 2 a (l — cos(t))g — esin(t)h, 


(5) 


e t adx e(ff) 


= ss 2 a (1 - cos(i))/ + (4 + c 2 a cos(t))g + an(t)h, 


(6) 


e* &dx ' (h) 


= es 2 a sin(i)/ - c 2 a sm{t)g + cos(t)h, 


(7) 




= (cos(2a) + 2s 2 cos(i))/ + e(cos(2o;) - 2c 2 cos(t))g - 2e sin(i)ft, 


(8) 


e tx e 


= 7rf + sin(t)X e + (l-cos(t))X 2 , 


(9) 


e tadf (g) 


sin(tc Q ) 
= cos{tc a )g H ft. 


(10) 



2.2 Existence of Optimal Trajectories 

A control system is complete if, for every measurable control function u : [a, b] — > [—1,1] and every initial 
state p, there exists a trajectory 7 corresponding to u, which is defined on the whole interval [a, b] and satisfies 
7(a) = p. Since 50(3) is compact and the function T{x, u) := x(f + ug) is regular enough, the system (1) is 
complete. Note that (f,g) satisfies the Strong Bracket Generating Condition (cf. [31]) and the set of velocities 
V(x) := {x(f + ug), u G [—1, 1]} is compact and convex. Then, (cf. for instance [36]): 

Proposition 1 For each pair of points p and q belonging to 50(3), there exists a time optimal trajectory 
joining p to q. 



2.3 Pontryagin Maximum Principle and Switching Functions 

We next state the Pontryagin Maximum Principle (PMP) (cf. [26]) for our minimum time problem on 50(3). 
Define the following maps called respectively Hamiltonian and minimized Hamiltonian : 

H : T*50(3) x [—1, 1] — > K, H(p,x,u) := < p,x(f + ug) >, (11) 
H : T*SO(i) — > K, H{p,x):= min H(p,x,v). (12) 

«€[-l,l] 

The PMP asserts that, if 7 : [a, b] — > 50(3) is a time optimal trajectory corresponding to a control u : [a, b] — > 
[—1,1], then there exists a nontrivial Geld of covectors along 7, that is an absolutely continuous function 
A : t G [a, b] X(t) G T^50(3) (identified with so(3)) never vanishing and a constant A > such that, for 
a.e. t G Dom{i), we have: 

A W = -f£(A(i),7(*)X*)) = -A(t)(/ + «(t)fl), 

ii) W(7(i),A(t),u(i)) + Ao = 0, 

iii) W( 7 (t),A(t),«(t))=H(7(t),A(t)). 
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Remark 1 The PMP is just a necessary condition for optimality. A trajectory 7 (resp. a couple (7, A)) 
satisfying the conditions given by the PMP is said to be an extremal (resp. an extremal pair). An extremal 
corresponding to Ao = is said to be an abnormal extremal , otherwise we call it a normal extremal . For a 
normal extremal, we can always normalize A = 1, and we do this all along the paper. Notice that in general an 
extremal corresponds to more than one covector. For this reason, usually, one distinguishes between abnormal 
extremal that arc strict (i.e. they correspond only to covectors satisfying Ao = 0) and abnormal extremal that 
are non-strict (i.e. they correspond to covectors with Ao = and to covectors with A ^ 0). 

A control u : [a, b] — > [—1, 1] is said to be bang-bang if u(t) e {—1, 1} a.e. in [a, b}. Moreover, if u(t) e { — 1, 1} 
and u{t) is constant for almost every t G [a, b], then u is called a bang control. A switching time of u is a time 
t e [a, b] such that, for each for every s > 0, u is not bang on (t — e, t + e) n [a, b] . A control with a finite number 
of switchings is called regular bang-bang. A trajectory of S is a bang trajectory, bang-bang trajectory, regular 
bang-bang trajectory respectively, if it corresponds to a bang control, bang-bang control, regular bang-bang 
control respectively. The switching functions, associated to an extremal pair (7, A), are the three "components" 
of the covector X(t) on the basis {/, g, h} transported to the point ^(t). More precisely: 

Definition 2 (switching functions) Let $i(x,p) (i = 1,2,3) be the Hamiltonian functions corresponding 
respectively to the vector fields f,g,h (cf. [18]). i.e. $i(x,p) :=<p,xf >, $>2{x,p) :~<p,xg >, $3(0;, p) .=< 
p, xh > and (7, A) be an extremal pair. The switching functions associated to (7, A) are the evaluations of 
&i(x,p) along the extremal i.e.: 

pi(t) := *i(7(*),A(t))=<A(t),7(t)/>, (13) 

■= *2(7(*),A(i))=<A(t),7(t)ff>, ( 14 ) 
¥>s(t) := $ 3 (7(*),A(t))=<A(t),7(t)/i>. (15) 

Remark 2 Notice that the <pi 's are at least continuous and since A never vanishes, the three switching functions 
cannot be all zero at the same time t. Moreover, using the switching functions, ii) of PMP reads: 

W(A(i),7(t),u(t)) = <pi(t) +u(t)ip 2 (t) + A = a.e. (16) 

The switching functions arc important because they determine where the controls may switch. In fact, using 
the PMP, one easily gets: 

Proposition 2 A necessary condition for a time t to be a switching is that (f2(t) = 0. Therefore, on any 
interval where if 2 has no zeroes (respectively finitely many zeroes), the corresponding control is bang (respectively 
bang-bang). In particular, Lp 2 > (resp ip 2 < 0) on [a,b] implies u = — 1 (resp. u = +1) a.e. on [a,b]. On the 
other hand, if f>2 has a zero at t and ^(t) exists and is different from zero, then t is an isolated switching. 

As a corollary, it holds a.e. along an extremal trajectory that: 

u(t)<p2{t) =-\M*)\- (1 7 ) 

An extremal trajectory 7 of £ defined on [c, d] is said to be singular if the switching function ip2 vanishes on 
[c, d]. To compute the control corresponding to a singular trajectory, one should compute the derivatives of the 
Pi's. Using the Lie bracket relations between f,g,h, one gets the system of differential equations (called the 
adjoint system) satisfied a.e.: 

if! = -uip 3 , (18) 
<f2 = ¥3, (19) 
lf 3 = sluif! - c 2 a tp 2 . (20) 

From Eqs. (15) and (19), one immediately gets that ip2 is at least a C 1 function. Moreover, if 7 is singular in 
[a, b], then ip 2 — and, from Eq. (19), we get ip 3 = a.e.. From (16) (cf. PMP ii)), we get that ipi = —1 a.e. 
on [a, b]. From (20) we get u = a.e. i.e.: 

Proposition 3 For the minimum time problem for (S), singular trajectories are integral curves of the drift, 
i.e. they correspond to a control a.e. vanishing. 
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In the sequel, we will use the following convention. The letter B refers to a bang trajectory and the letter S 
refers to a singular extremal trajectory. A concatenation of bang and singular trajectories will be labeled by 
the corresponding letter sequence, written in order from left to right. Sometimes, we will use a subscript to 
indicate the time duration of a trajectory so that we use B t to refer to a bang trajectory defined on an interval 
of length t and, similarly, St for a singular trajectory defined on an interval of length t. 

If we fix u € [—1,1], then the integral curves of x(f + ug) are periodic. In particular, the integral curves of 
xX e are periodic with period 2tt while the integral curves of the drift xf are periodic with period 2it/c a . This 
means: 

Proposition 4 If 7 is an extremal trajectory of type B t (resp. S t ), then t < 2n (resp t < 2^jc a ). 

There are two quantities that are remain constant along an extremal trajectory. The first one comes from the 
fact that the minimized Hamiltonian H is constant along the extremal pairs (7, A) (cf. (16) and 17): 

h := -¥>i(t) + \<P2(t)\ = A o, (21) 

with A equal to zero or one (cf. Remark 1). The second conserved quantity is: 

h := c 2 a <pl + s 2 a tpl +<pj = K 2 , for some K e K. (22) 

Remark 3 Equations (18), (19), (20) are Hamiltonian equations on the dual of so(3), with respect to the 
canonical Poisson structure induced by the brackets of /, g, h e so(3), and corresponding to the left-invariant 
Hamiltonian (11). The conserved quantity I2 is the Casimir function (see for instance [1]). 

There is a geometric interpretation of the above equations. Let (7, A) be a normal extremal lift of the 
time-optimal control problem. Then, the adjoint vector A with coordinates (<£i)i=i,2,3, lies in the intersection 
of the region defined by Eq. (21) and the ellipsoid defined by Eq. (22). 

2.4 Classification of optimal trajectories 

In this section, we investigate the structure of time optimal trajectories by analyzing the extremal flow defined 
in (18)-(20), subject to (21) and (22). First we study abnormal extremals (we prove that they are regular 
bang-bang and we establish a relation between the interior bang times) . Then we study normal extremals that 
are bang bang (again we find relation between the interior bang times). Finally we study optimal trajectories 
containing a singular arc. The results presented in this section are well-known, and some of them already 
contained in [3, 4], although in many cases without proof. To have a self-contained paper, we provide an 
argument for all of them. 

2.4.1 Abnormal Extremals 

The following proposition describes the switching behavior of abnormal extremals. 

Proposition 5 Let 7 be an abnormal extremal. Then, it is regular bang-bang and the time duration between 
two consecutive switchings is always equal to ir. In other words, 7 is of kind B^B^.^B^Bt with t <ir. 

Proof of Proposition 5 By definition, Ao = 0. Then Eq. 16 becomes 

tpi(t) = -u(t)<p2{t), for a.c. t € Dom(j). (23) 

If 7 is singular on some interval [c,d], then ip 2 = and from (19) (^3 = on [c, d]. Eq. (23) gives ipi = 0, 
contradicting the non triviality of A (cf. Remark 2). Then 7 cannot contain a singular arc. Therefore, u 2 = 1 
a.e. t G Dom("f). 

From Eqs. (20) and (23), we get a.e. (ps(t) = {—s 2 a u(t) 2 — c 2 a )ip2{t) = —<p2{t). This means that, in the 
(ips, (p 2 ) plane, the vector z{t) :— (<fi3(t), y>i{t)) rotates with angular velocity equal to one (cf. Eq. (19)). This 
implies 7 is a regular bang-bang trajectory and the time duration between two consecutive switchings along 7 
is always equal to ir. H 
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2.4.2 Normal Bang-Bang Extremals 

Let 7 be a bang-bang trajectory starting at po and ending at p e^ taX ^ e^ tlX -^ e^ t2X +^ e^ tsX -\ The case in which 
the first bang is of kind A_ is similar. We have 932(^0) = ^2(^0 + ti) = ^2(^0 + t\ + t 2 ) = which implies: 

< \{t + t 1 ),p 2 e' t ^ adx -(g) >=< X(t +t 1 ),p 2 g >=< X(t + h) 1 p 2 e t - adx + (g) >= 0, (24) 

where p 2 = pae^ toX+ ^e^ tlX -\ We need the following definition. If z\, z 2 , z 3 are (possibly time-varying) vector 
fields of 50(3), the application q 1— > q{z\ A z 2 A 2:3) is the field of 3-vectors associated to the z^s, where 
q(zi A z 2 A z 3 ) is an element of /\ T q SO(3), the 3-fold exterior power of T q SO(3). We now rewrite Eq. (24) by 
using fields of 3- vectors. We obtain: 

g A e ~ tiadx - (g) A e t2adx + (g) = 0. (25) 

Thanks to (6), Eq. (25) is equivalent to r(t\,t 2 )f A g A h = for an appropriate real-valued function r. After 
computations, we get r{t\,t 2 ) = sin( * 1 ~* 2 ). This implies that t\=t 2 = t 3 . 

Similarly to what we did in the proof of Proposition 5, consider now a time optimal trajectory of the form 
BBtB, where Bj- is a nontrivial interior bang arc associated to a normal extremal and T g]0, 2ir[. From (20) 
and (16) (with Ao = 1), we get 933 = —(<p 2 + s^u). Using (19), this means that the vector z = (953, ip 2 + us"^) T 
satisfies the differential equation 




z, te]o,T[, 



with boundary conditions (the switching conditions imply (p 2 (0) = (fi 2 (T) = 0) z(0) = (933(0), ms^) and z(T) = 
(ifi3(T),us^). Using 932(0) = and the fact that u > (resp. u < 0) implies > ^2(0) = <p 3 (0) (resp. 
< 932(0) = 933(0)), one easily gets tan(T/2) = 9? 3 (0)/(u s 2 a ) < 0. It follows that T e {tt,2tt). In summary, we 
proved that 

Proposition 6 Let 7 be a bang-bang normal extremal. Then the time duration T along an interior bang arc 
is the same for all interior bang arcs and verifies tt < T < 2tt. 

Remark 4 From Propositions 5 and 6, we get that, for an extremal bang-bang trajectory (normal or abnormal), 
the time duration T along an interior bang arc is the same for all interior bang arcs and verifies tt < T < 2tt. 

2.4.3 Optimal Trajectories Containing a Singular Arc 

The purpose of this paragraph is to describe the structure of time optimal trajectories containing singular arcs. 

Proposition 7 Let 7 be a time optimal trajectory containing a singular arc. Then 7 is of the type B t S s B t ' , 
with s < if t > or t' > and s < 2^- otherwise. 

Proof of Proposition 7 Let 7 be a time optimal trajectory containing a singular arc St, t > 0. From 
Proposition 3, we know that t < 2ir j c a . 

Assume now that 7 contains a singular arc and a nontrivial interior bang arc. Then, we may assume that 
7 contains a piece of the type S s B t or B t S s (say the first), with B t a complete bang arc. Then we have 
^2(5) = 93 3 (s) = 9? 2 (s + t) = 0. This translates to: g A h A e tadx '~ (g) = 0. Using (6), it implies that cos(f) = 1, 
i.e. t = 2ir. This contradicts the time optimality of 7. Finally, from (10), we get e~ (g) = —g. From this, we 
deduce for t > 0, e^ f e tXe — e tx - E e^ f . Therefore, for t, s > 0, we have e sf e tXc = eS s ~^^ e tx - e e^^ . Then, 
if s > 2- and t > and taking into account what precedes, e s) 'e tXe cannot be optimal. H 

2.5 Uniform Bound on the Number of Switchings for Time Optimal Trajectories 

For a e]0, 7r/2[, let N(a) be the supremum of the number of switchings of any time optimal trajectory on 
SO(3). Thanks to the left-invariancc of the control system (1), we may assume that the supremum is taken 
over any time optimal trajectory starting at Ld. In this paragraph, we prove the following: 

Proposition 8 For a e]0, ir/2[, N(a) is finite (and thus achieved). 
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Proof of Proposition 8 Let us first prove that: 

Claim every optimal trajectory of (£) is a finite concatenation of bang and singular arcs. 

Let 7 : [a, b] — > 50(3) be a time optimal trajectory of (£). Let 5 be the set of zeroes of ip2 such that, if t € 5, 
then ip 2 does not vanish identically in some neighborhood of t. Clearly 5 is the set of times t such that j(t) 
is the junction of two bang arcs or the the junction of a singular arc and a bang arc. The conclusion follows 
if 5 is finite. Reasoning by contradiction, 5 must have a limit point t. Moreover t € 5, otherwise (f2 would 
vanish identically in a neighborhood of i, contradicting the fact that t is a limit point of 5. Note also that 
ip 2 is continuous in an open (in [a, b]) neighborhood N of t (see Remark 2). By definition of 5 and t, there 
exists a sequence (i„) in N converging to t such that (p 2 {t n ) 7*= 0. Pick n large enough, so that, if [t' n , t'£\ is the 
maximal subinterval containing t n with tp 2 7*= on (t' n , i"), then [t' n , t n ] C N. Clearly (p 2 (t' n ) = ¥>2{t n ) = 0, 7 is 
a bang arc on [t' n , t n ] and — t' n tends to zero as n goes to infinity since t n tends to F. But, by Proposition 6, 
t'n — t' n > 7T for n large enough. So we reached a contradiction and 5 is finite. The claim is proved. □ 

To finish the proof of Proposition 8, it remains to show that the (finite) number of switchings for any 
time optimal trajectory is uniformly bounded over 5*0(3). The argument goes by contradiction: there would 
exist, then, a sequence of regular bang-bang time optimal trajectories B Sn Bx n • • ■ BT n B tn , where s n ,t n < 2n, 
7r < T n < 2ir and the number of switchings m„ goes to infinity as n goes to infinity. Therefore, there exists a 
sequence of points (x n ) of 5*0(3) such that, the minimum time r„ needed to connect Id to x n by a trajectory 
of (E), goes to infinity as n goes to infinity. 

To reach a contradiction, it is enough to show that there exists a time T so that, for every point x G 50(3), 
there exists a trajectory 7 of (£) connecting Id to x with T(j) <T. By a compactness argument and thanks 
to the 50(3)-invariance of (E), that would result from the following fact: there exists t > and an open 
neighborhood U C 50(3) of Id such that every point x € U can be reached from Id in time less or equal than 
t. The latter simply results from the facts that (S) has the accessibility property and e 1 ' is periodic. H 

Remark 5 Since the degree of non-holonomy of the distribution generated by (/, g) is equal to two, then, by 
standard controllability arguments, one can quantitatively relate the size of U and t as follows: U contains a 
ball of radius S| for some positive constant C. Therefore, N(a) can be bounded above by for some positive 
constant C. 



3 The Hopf Fibration and Proof of Theorem 1 
3.1 The Hopf Projection 

In this paragraph, we describe explicitly the Hopf projection from 50(3) to MP 2 . This projection provides 
50(3) with a structure of fiber bundle with base RP 2 and fiber 5 1 . In the sequel, we use the identification of 
KP 2 with 5 2 \ <~, where ~ is the antipodal map, that is RP 2 is the set of rows (2/1,1/2,2/3), J2vf = ^ where 
(2/1,2/2,2/3) ~ ( — j/i, — 2/2, — 2/3)- I n the sequel, RP 2 is identified with the subset NI , made of the (open) top 
hemisphere together with half of the equator. Fix a point yo G RP 2 . The Hopf projection is defined as: 

n : 50(3) -» RP 2 , 

x >-> 2/ = Vox, 

where yox is the standard matrix product. Then, any left-invariant vector field V : x xv on 50(3), v € so(3), 
is transformed by dH into the (left-equivariant) vector field Vs = dIL(V) : y yv. In the following, we call 
respectively, the control systems x = x(f + ug), x e 50(3) and y = y(f + ug), y e RP 2 , the control systems 
upstairs and downstairs . As explained next, it is crucial to choose yo so that the drift term at the initial point, 
Id /, vanishes downstairs, so we require that IL(Id) = yo- Indeed, if it is the case, then: i) from the point yo, we 
have local controllability (see next section) and this greatly helps in the construction of the optimal synthesis; 
ii) if a trajectory upstairs starts with a singular arc (that is, with u = 0, i.e. it is an integral curve of the drift 
/), then its projection is a point. That suitable choice of yo is made possible by the following normalizations: 





1 


o\ 


1 { 















and g — s a 












J 


V 


-1 


/ 



f = ca -1 and<7 = s Q 1 , y = (0,0,1). 
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In that way, the control system downstairs reads: y = Fs(y) + uGs(y), where Fs(y) — yf, Gs(y) — yg. In 
order to respect the convention in control theory where states are represented by column vectors and costates 
by row vectors, we will consider the transposed control system and, with the change of notations y T — > y, 
f T = — / — ► /, g T = —g — > g, we obtain downstairs: 



V = F s (y)+uG s (y), \u\ < 1, 

0-10 
1 


0-1 



Fs{y) = fy = fxy = c, 
G s {y) 



gy = gxy 



o l o 





(26) 



3.2 Proof of Theorem 1 



For the rest of the paper, we assume a g]0, 7r/4[. In this Section we prove Theorem 1, using a Lemma describing 
the structure of time optimal trajectories of (£)s connecting the north pole to any point of Nl, that will be 
studied in the next section. 

Lemma 1 Consider the control system (S)s and the time optimal trajectories connecting the north pole to any 
point of Tff. Then 

i) they are regular bang-bang with same time durations for the interior bang arcs, i.e., they are of the type 



B X B, 



v(s) 



B v ( s )B t , where s e [0,7r] ; t G [0, v(s)] and 

v(s) = tt + 2 arctan(- 



cot 2 (a) 

ii) Set N (a) := 2 [^] . Then the maximum number of switching of these trajectories is 



N s (a) = N (a) - JV (a) 



4a 



(27) 



(28) 



The above formula means that JVs(a) can take the values N (a), N (a) + 1 or N (a) + 2; 
iii) They are projections, through the Hopf map H, of time optimal trajectories of (£) starting at Id. 
Proof of Lemma 1 

For the Proof of i), see Proposition 10 and Proposition 11. For the proof of iii), see Lemma 7 and Lemma 8. 
For the proof of ii) see Proposition 13. H 



We now prove separately the two inequalities of Theorem 1. 

Proof of the inequality Ns(a) < N(a). From iii) of Lemma 1, every time optimal trajectory of (£)s 
connecting the north pole to any point of N£ is the projection by II of a time optimal trajectory of (£) with 
the same time duration (in particular of the time optimal trajectory connecting the two fibers). Therefore, 
N s (a) < N(a). □ 

Proof of the inequality N(a) < Ns(a) + 4. We refer to Fig.l. Consider a time optimal trajectory 7 of (S) 
containing N(a) switchings. With no loss of generality, we may assume that N(a) > 2. By Propositions 6 and 
7, we deduce that 7 is regular bang-bang and is of the type B s Bt • • • BrB t , with s, t > 0, tt < T < 2n. Since 
every subarc of a time optimal trajectory is also time optimal, we may assume that s = t = 0. Let Id and x\ 
be the initial and terminal points of 7 and consider 75, a time optimal trajectory for (£)s connecting U(Id) 
and II(xi). From i) of Lemma 1, 75 is of the type B s /B v t s ^ ■ ■ ■ B v ( s ,}B t > with s' < tt, t' < v(s') and m interior 
bangs. We thus have m < Ns(a) — 1. We now build, from 75, a suboptimal trajectory connecting Id and x\ 
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Figure 1: Proof of Theorem 1 



as follows: we can lift 75 to SO(3) to an admissible trajectory 75 of (S) connecting 20 and x\, with xo in the 
fiber of H(Id). It is also clear that 7s and js have same time durations. By construction of the fiber of H(Id), 
we get that xo = e* ' with t" < 2n. Finally, the curve 7 obtained as the concatenation of e* ^ and 75 is an 
admissible trajectory of (£) connecting Id and xi. Its time duration is equal to 

T(7) = i" + T(7s) = i" + T( 7s ) = t" + mv(s') + s' + t', 

with m < Ns{a) — 1. Since 7 is time optimal, we have T("f) < Ttf), which implies that 

(N(a) - l)T < t" + mv(s') + s' + t'. 

Using all the estimates on T,s',t',t" (i.e. T e [w, 2ir), cf. Remark 4, s' < w, t' < max s / e [(|.j] v(s'), t" < 2ir), 
we deduce that 

( N(a) — 1 )ir < Ns(a)V(a) + 37r, where V(a) := max se [ 0;7r ] v(s), 

from which we have 

N(a) - N s (a) < N s (a) V ^ - - + 4. (29) 

Set r(a) := iVg(a) v ^ t ^~ ,r . A simple computation shows that r(a) = Ns(a)^ arcsin(tan 2 (a)). Using (28), it is 
easy to see that r(a) e]0, 1[ on ]0, 7r/4[. Since iV(a) and Ns(ce) are integers, we get N(a) — Ns(a) < 4. □ 



4 The Time Optimal Synthesis Downstairs 

In this section, to compute Ns(a), we study the time optimal synthesis for the problem downstairs (26), starting 
from the point yo- 

Definition 3 A time optimal synthesis for the problem downstairs (26), starting from the point yo is a family 
of time optimal trajectories T — {j y : [0, b y ] \— ► RP 2 , y e RP 2 : 7 y (0) = y , "f y {b y ) = y}. 

For that purpose, we use the theory of optimal syntheses on 2-D manifolds developed by Sussmann, Bressan, 
Piccoli and the first author in [9, 10, 11, 14, 15, 23, 24, 32, 33] and recently rewritten in [12]. The core of 
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the Theory consist of an explicit algorithmic construction (by induction on the number of switchings) of the 
optimal synthesis. 

Note that the previous theory uses a a more elaborated concept of synthesis, namely that of regular synthesis 
(see for instance [8, 12, 16, 25] and cf. Section 4.1.1). In the following, in order to compute Ng(a) we just 
need to follows the steps of the algorithmic construction mentioned above, without requiring the existence of a 
regular synthesis. In the sequel, by time optimal synthesis, we refer to one in the sense of Definition 3, whose 
existence is simply guaranteed by Proposition 1. 

Consider a two dimensional smooth manifold M and the problem of computing the time optimal synthesis 
from a fixed point yo € M for the control system: 

y = F(y) + uG(y), y e M, \u\ < 1, (30) 

where F and G are C°° vector fields. We introduce three functions: 

A A (y) := Det(F(y),G(y))=F 1 (y)G 2 (y)-F 2 (y)G 1 (y), (31) 

A B (y) := Det{G{y) 1 [F,G]{y))^G l {y)[F,G] 2 {y)-G 2 {y)[F,G] l {y), (32) 
f s (y) := -A B (y)/A A (y). (33) 

The sets A^(0), A^O) of zeroes of A a, Ab are respectively the set of points where F and G are parallel, and 
the set of points where G is parallel to [F, G] . These loci are fundamental in the construction of the optimal 
synthesis. In fact, assuming that they are smooth embedded one dimensional submanifold of M we have the 
following: 

• in each connected region of M \ (A^ : (0) U A^ 1 (0)), every extremal trajectory is bang-bang with at most 
one switching. Moreover, if the trajectory is switching, then the the value of the control switches from 
-1 to +1 if f s > and from +1 to -1 if f s < 0; 

• the support of singular trajectories (that are trajectories for which the switching function identically 
vanishes, see Definition 4 below) is always contained in the set A^ 1 (0); 

• a trajectory not switching on the set of zeroes of G is an abnormal extremal (i.e. a trajectory with 
vanishing Hamiltonian) if and only if it switches on the locus A^ 1 (0). 

Then the synthesis is built recursively on the number of switchings of extremal trajectories, canceling at each 
step the non optimal trajectories (see [12], Chapter 1). 

Remark 6 As we will see later (see Proposition 10), the condition that a < tt/A guarantees that there are no 
singular trajectories for the problem downstairs. 

4.1 Basic Definitions and Facts on Optimal Synthesis on 2-d Manifolds 

Consider the minimum time problem for the control system (30). In this Section, we recall some key facts for 
the construction of time optimal synthesis following [12]. 

The first ingredient is, as usual, the PMP, that, on a two dimensional manifold, has exactly the same form as 
described in Section 2.3 but with the following change of notation: x € 50(3) — ► y E M, A(t) € T^qSOfi) — > 
X(t) € T 7 ( t )M. As for the problem upstairs, switchings are described by the switching function: 

Definition 4 (Switching Function) Let (7, A) be an extremal pair. The corresponding switching function is 
defined as (f>(t) :=< A(t), G(-y(t)) >. 

Again, 4> is at least continuously differentiable (4>(t) =< \(t), [F, G] {"f(t)) >, cf. discussion in (19), and it 
determines the switching rule, according to Proposition 2 with the change of notation ip 2 — > <f). Again, an 
extremal trajectory 7, defined on [a, b], is called singular if <j> = in [a, b]. The following three Lemmas 
illustrate the role of the two functions defined in (31), (32). The proofs can be found in [9, 12, 24]. 
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Figure 2: Lemma 5 and Lemma 6 



Lemma 2 Let 7 fee an extremal trajectory that is singular in [a, b] C Dom(-f). Then 7|[ a ,b] « s associated to the 
so called singular control <p{"i{t)), where: 



viy) = 



VA B (y)-F(y) 
VA B (y)-G(y)- 



(34) 



with Aa and Ag defined in Eqs. (31) and (32). Moreover, on Supp(j), ip(y) is always well-defined and its 
absolute value is less than or equal to one. Finally Supp(-j\[ a .b]) C A iJ 1 (0). 

Lemma 3 Let 7 be an extremal bang-bang trajectory for the control problem (30), to G Dom(^) be a time such 
that 4>{to) — and G(j(to)) 7^ 0. Then, the following conditions are equivalent: i) 7 is an abnormal extremal; 
ii) j(t ) e A^O); iii) ~/{t) e A^O), for every time t e Dom{^) such that <f>(t) = 0. 

The following lemma describes what happens when and As are different from zero. 

Lemma 4 Let fl C M be an open set such that Q n (A^ 1 (0) U Ag 1 (0)) = 0. Then all connected components of 
Supp{^)C\£l, where 7 is an extremal trajectory of (30), are bang-bang with at most one switching. Moreover, if 
fs > throughout Q, then j\q is associated to a constant control equal to +1 or — 1 or has a switching from 
— 1 to +1. If fs < throughout 0, i/ien 7|n is associated to a constant control equal to +1 or —1 or has a 
switching from +1 to —1. 



Definition 5 Let 7+ : [0,r] — > M (resp. 7~ : [0,r] - 
corresponding to the constant control u = 1 (resp. u = 
support of the curve 7 + |[o.t] fresp. 7[Q t jj- 



M) be the trajectory of (30) starting at y and 
-I). Fort e]0,r] ; Zet fresp. r~(i); 6e tfie 



Under the assumption F(y ) = and Ag(j/o) 7^ 0, the next lemma (for a proof, see for instance [12, 26]), 
describes the shape of the optimal synthesis in a neighborhood of y . That local behavior of the optimal 
synthesis remains actually the same as long as T + (t) and T~(t) do not intersect A^ 1 (0) and A^ 1 (0) (except 
of course at y ). 



Lemma 5 Consider the control system (30). Assume that F(y ) = and As(yo) 7^ 0. Let £1 be an open 
neighborhood of yo such that, Q n A i3 1 (0) = and O n A^ : (0) is an embedded one- dimensional submanifold of 
Q. Let 7+ : [0,t] -> M (resp. 7" : [0,r] -> MJ 6e ifte trajectory of (30) start ing at yo and corresponding to 
the constant control u = 1 (resp. u = — lj. Then, for every t+,t- g]0,t[ smc/i t/ia£; (a) r + (t + ), T~(i_) C O, 
(b) T+(t+) n A A : (0) = r-(t_) n A^(0) = {y }, (c) r+(t+) n r~(t_) = {y }, we ftaue the following. There 
exists an open neighborhood U ofT + (t + ) U T~(t-) contained in £1 such that, for every y G U, there exists a 
unique extremal trajectory of (30) of the type B s B t contained in U , which is time optimal and steers yo to y. 
In particular, the system (30) is controllable in U and 7+ (resp. 7~ ) is time optimal up to t + (resp. t-), see 
Fig. 2 A. 
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Finally, we need one more lemma, related to Lemma 3, and whose hypothesis are illustrated in Fig. 2 B: 

Lemma 6 Consider the control system (30). Assume that i) F(y ) — 0, Ab{vo) 0; ") there exists t + > 
such that r+(t+) n A A X (0) = {y , 7+ iii) there exists e > such that T+(t + + e) n A^O) = 0. Then 

7+ is extremal exactly up to time i + . Moreover, any extremal trajectory 7 defined on [0, T] with T > t + and 
coinciding with 7+ on [0, t + ], switches at t + to the constant control u = — 1 one? i/ius 7 is an abnormal extremal 
(cf. Lemma 3). A similar statement holds for j~ . 

Remark 7 Under the hypotheses of Lemma 6, one can prove that the abnormal extremal 7 restricted to an 
interval [0,T], is a non-strict abnormal extremal if T < i + , while it becomes a strict abnormal extremal if 
T >t + (cf. Section 2.3). In other words, 7 becomes a strict abnormal extremal after the first switching. These 
fact are analyzed in details in [11] and [12] (see Chapter 4, and in particular Section 4.3, where strict abnormal 
extremals are called Non Trivial Abnormal Extremals). 

4.1.1 Frame Curves and Frame Points 

In this paragraph, we briefly recall, for sake of completeness, the main results of the theory developed in 
[23, 24] (see also [12]). That material is only used here and in Section 5, where some numerical simulations and 
conjectures are presented. In [23, 24] (see also [12]), it was proved that the control system (30), under generic 
conditions on F and G (with the additional assumption F(y ) = 0) admits a time optimal regular synthesis in 
finite time T , starting from j/o- By generic conditions, we mean conditions verified on an open and dense subset 
of the set of C°° vector fields endowed with the C 3 topology (see [12], formula 2.6 pp. 39). More precisely, let 
1Z(T) be the reachable set in time T > given by: 

K(T) := {y e M : 3 b y e [0,T] and a trajectory 

7 y : [0, by] -> M of (30) such that j y (0) = y , %(b y ) = y}. 

Then a time optimal regular synthesis is defined by: i) a family of time optimal trajectories T = {j y : [0, b y ] — > 
M, y e TZ(T) : 7 y (0) = yo, ~f y {b y ) = y} such that if y y e T and y = j y ^ for some t e [0, b y ], then j y = 7 w |[ ,t]; 
ii) a stratification of 1Z(T) (roughly speaking a partition of TZ(T) in manifolds of different dimensions, see [12], 
Definition 27, p. 56) such that the optimal trajectories of T can be obtained from a feedback u{y) satisfying: 

• on strata of dimension 2, u(y) = ±1, 

• on strata of dimension 1, called frame curves (FC for short), u(y) = ±1 or u{y) = ip(y), where ip(y) is 
defined by (34). 

The strata of dimension are called frame points (FP). Every FP is an intersection of two FCs. In [24] (see 
also [12]), it is provided a complete classification of all types of FPs and FCs, under generic conditions. All the 
possible FCs are: 

• FCs of kind Y (resp. X), corresponding to subsets of the trajectories 7+ (resp. 7") defined as the 
trajectory exiting yo with constant control +1 (resp. constant control —1); 

• FCs of kind C, called switching curves, i.e. curves made of switching points; 

• FCs of kind S, i.e. singular trajectories; 

• FCs of kind K , called overlaps and reached optimally by two trajectories coming from different directions; 

• FCs which are arcs of optimal trajectories starting at FPs. These trajectories "transport" special infor- 
mation. 

The FCs of kind Y, C, S, K are depicted in Fig. 3. There are eighteen topological equivalence classes of FPs. 
A detailed description can be found in [10, 12, 24]. 
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Figure 3: 



Remark 8 The proof of the existence of a regular synthesis is made by means of a constructive algorithm 
(working recursively on the number of switchings), that builds explicitly the optimal trajectories (see [12], 
Section 2.5 p. 56). We stress the fact that the existence of a regular synthesis cannot be guaranteed before the 
complete execution of the algorithm. Since for our systm (26), we do not reach the end of that construction, 
we cannot conclude that such a regular synthesis exists. However, we conjecture that last fact (see also Section 
5). 

4.2 The Problem Downstairs 

In this section, we apply the theory recalled in Section 4.1 to the control system (26) on S 2 in order to compute 
Ns(ct), the maximum number of switchings for time optimal trajectories connecting the north pole to any point 
of jW. First we need some notations. 

Definition 6 Set: 

-c a y 2 

X+(y) = F s (y) +G s (y) = X + y = X + x y = | c a -s a ] ( y 2 ) = ( c aVl - s a y 3 

SaV2 

-c a y 2 

x s(v) = F s(y) ~ G s {y) = X~y = X~ x y = \ c a s a ) | y 2 | = | c a y x + s a y 3 

-s a y 2 

Let 7 : [t\,t 2 ] — > S 2 be a trajectory of (30). If 7 corresponds to the constant control +1 (resp. -1) in \t\,t 2 ], 
we say that ^\\ tl ^ 2 ] is a X + -trajectory (resp. X~ -trajectory). Moreover, we call 7 ± the trajectories exiting the 
point xo with respectively, constant control +1 and — 1. Let t^ be the last times for which 7 ± are optimal. We 
define 7^ := 7 ± |[ i±]- If Ji '■ [ a i b] — ► S 2 and 72 : [b, c] — > S are trajectories of (26) such that 71(6) = 72(b), 
then the concatenation 72 * 71 is the trajectory: 
















s 





•5 a 








c a 










^ a 

















) 




-[ 









(72 *7i)(*) : = 



71 (t) /ori g [a, b], 

72 (t) /orte [6,c]. 



Notice that, in the notation 72 * 71 , 71 comes first. 

The first quantities to be computed are A^ 1 (0), A^ 1 (0) and the sign of fs- Referring to Fig. 5, we have for the 
system (26): 

A^(0) = {( yil y 2 ,y 3 ) T e S 2 : y 2 - 0}, 
A^iO) = {( yil y 2 ,y 3 ) T e S 2 : y 3 = 0}, 

fs(y) > 0, V y g {(j/i,y 2 ,y 3 ) T e 5* 2 : y 2 y 3 > 0}, 

/ s (y) < 0, V y G {(2/1, ?/2, 2/3 ) T G 5 12 : 2/22/3 < 0}. 

(35) 
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Figure 4: 



The set A B X (0) is called the equator and A A 1 (0) the meridian. Moreover, let NI be the (open) top hemisphere, 
i.e. the set of points (yi,y 2 ,y3) T so that 2/3 > and (see Fig. 4): 

NI+ := {y e NI : y 2 < 0}, M+ = {y e Aff : Wl > 0, y 2 = 0}, £ + = {y£S 2 :y 2 < 0, y 3 - 0}. 

Similarly 

m~ = {y e NI : y 2 > 0}, M" = {y e : j/i < 0, y 2 = 0}, r = {ye5 2 : te > 0, y 3 = 0}. 

We also parametrize points y of the meridian by the oriented angle between Ot/o and Oy. We use P(£), 
£ G [— 7r,7r], to denote the point of the meridian defined by the angle £. Then -P(O) = yo and P(a) (resp. 
P(—a)) is the center of rotation in the north hemisphere of Xg (resp. Xg). We also have that j + (resp. 7"), 
up to time n, is a half-circle with diameter [y ,P(2a)] (resp. [j/0j -P( — 2a)]), see Fig. 4. 
From Lemma 4, it follows: 

Proposition 9 Let 7 : [0, T] — > S* 2 , 7(0) = yo &e an optimal trajectory for the control system (26). Then: 

• 7 /ias af most a * X - switching in NI~ , that is, if Supp(j\[ a ^) C NI~ , then j\[ a ^] corresponds to 
one of the three following controls: 

(-) u = +1 in [a, b], 
(-) u = — 1 in [a, b], 

(-) i/iere exists c <s]a, 6[, smc/i i/iai u = — 1 in [a, c[ and u = +1 in ]c, 6]; 

• 7 has at most an Xg * Xg switching in Nt + ; 

• 7 has at most an Xg * Xg switching in the region {x G S 2 : y 2 > 0, j/3 < 0}; 

• 7 has at most an Xg * Xg switching in the region {x € S 2 : y 2 < 0, j/3 < 0}. 

In Fig. 5, the integral curves of Fs, Gs, Xg , Xg and the loci A^ 1 (0), A^ 1 (0) are depicted. Moreover, the 
allowed switchings are indicated. 

Remark 9 Notice that, in Nt + (resp. NJ~), Xg points on the right (resp. on the left) of Xg , while, on the 
meridian, Xg and Xg are parallel (see Fig. 4). More precisely, Xg and Xg point in the same direction on 
{-P(£)> £ S]a, 7T— a[ U ] — 7T + o:, — a[} and in opposite directions on {P(£), £ €] — a, a[ U ]n — a, ir] [J [—it, —n + 
«[}■ 
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4.2.1 Two properties of extremal trajectories 

The following two propositions are essential in the construction of the optimal synthesis. 

Proposition 10 Every time optimal trajectory of (26), starting at the north pole, is regular bang-bang. 

Proof of Proposition 10 Since a < j,by taking into account Lemmas 5 and 6, the curves 7 + and 7~ defined 
in Definition 6, do not intersect the equator and are time optimal until the first time they meet the meridian, 
i.e. exactly up to time tt. Moreover, since singular arcs are contained in the equator and, thanks to Lemma 6, 
any time optimal trajectory 7 of (26), with at least one switching, is of the form B s B t ..., with s G]0,7r] and 
t > 0. Finally, since 7 is the projection of a time optimal trajectory 7 of (1), then the latter is also of the type 
B s B t .... Therefore, by Proposition 7, 7 cannot contain any singular arc, and so 7. H 

Proposition 11 Let 7 : [0,T] — > S 2 be a time optimal trajectory for the control system (26) of the type 
B s Bt 1 B t2 .... Then, all time durations of interior bang arcs are equal to v(s), where: 

v(s) := 7r + 2 arctan ( Ss „, - ) . (36) 

\c s +cot 2 (a)J 

Proof of Proposition 11 Consider 7 : [0,T] — > 5*0(3), an optimal trajectory that projects on 7 through the 
Hopf fibration EL Thanks to Proposition 6 (see also Remark 4), we have 7 = B s B tl B tl • • •, where t\ £ [tt, 2tt[. 
Moreover, since that curve projects on a time optimal trajectory for (26), we will establish a relation between 
s and t\. We start from the relations <pz{s) — <£i{s + t\) = 0, which can be written 

< X(s), Gs(j(s)) >=< A(s + h), Gs(j(s + ii)) >= 0. (37) 

Recall that A(s) = \{0)e- sX ' , X(s + h) = \(0)e- sX 'e~ tlX - and 7(a) = e sX ^{0), j(s + ti) = e tlX -- e sX ^(0). 
Since 7 is nontrivial, then A(0) is a nonzero line vector of R 3 . Moreover, 7(0) = yo = (0, 0, 1) T . Eq. (37) can 
be written as: 

X(0)e- sX '(g x e sX - 7 (0)) = 0, A^e^e"* 1 *-^ x e* lX -e sX " 7 (0)) = 0. 
The previous equations can be transformed to: 

det(e s ^A(0) T ,.g,e s ^ 7 (0)) = 0, det(e tlX -^ e sX ^ A(0) T , g, e tlX -^ e sX ^(0)) = 0, 

and then to: 

det{e sX ^\(0) T ,g,e sX 'j(0)) = 0, det(e sX ^ A(0) T , e~ tlX ^g, e sX ^(0)) = 0. 
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Since e sXc X(0) T is not zero, we deduce that: 

dct( 3 ,e sJ S(0), e - tlX -^) = 0. (38) 

We end up with the relation: 

-si cos(s - h/2) = c 2 a cos(ii/2). (39) 
Taking into account that ir < t\ < 2ir, we can simplify the previous equation to get (36). H 



4.3 Construction of the time optimal synthesis 

In this section, we present, step by step, the construction of the TOS for (26). Since we will not complete that 
construction, we only provide here the steps for which the outcome is justified by a rigorous argument. For the 
other steps of the construction, we refer to the last section where we propose conjectures on their outcomes, 
which arc supported by numerical simulation. 

Step 1 By Lemmas 5 and 6, for every e > 0, there exists an open neighborhood U of r + (7r — e) U r~(7r — e) 
(recall Definition 5) where the time optimal synthesis is described in Fig. 6 A. Moreover, £+ p = t~ p = tt 
(recall Definition 6); 

Step 2 taking into account the analysis of Sections 4.1 and 4.2, the time optimal trajectories for the problem 
downstairs are described by the following: 



Proposition 12 Every time optimal trajectory for the system (26), starting from the north pole, is 
contained in the following two sets of extremals, which are parametrized by the length of the first bang 
arc, the one of the last bang arc and the number of arcs: 

m terms 

E+(s,t) = ^V^W.-.^'We^'jo, (40) 
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S-(s,t) = ef£VVj^ ■ ■ ■ eff^Wef^ j/o, (41) 
m' ierms 

w/iere s G [0,7r], i e [0,u(s)] ; the number of bang arcs (m and m! respectively) is an integer and 

(-) e = +1 (resp. £ = —I), if m is odd (resp. even), 
(-) e' = +1 fresp. e' = —I), if m' is even (resp. odd). 

Step 3 Let A + and A~ be the two extremal trajectories starting resp. with controls u = 1 and u = —1, and 
switching after time n, i.e. corresponding resp. to E + (tt,-) and S~(7r, •). These two curves are abnormal 
extremals and their respective first bang arcs coincide with 7+, and ^y~ p . As explained in Remark 7, these 
two curves become strict abnormal extremals after time it. 

To describe them, consider, for e = ± and < k < k (k defined below), the half-circles L| C Clos(NI e ), 
whose centers lie on OP(ea) and passing through the points P k e and P§ +1 , where 

P+ := P(2na), P~ := P(-2na), 

for the integers n so that 2na < § + 2a. Note that ^ + 2a < it for a < j and, in fact, the last P r > belongs 
to the bottom-half hemisphere, i.e. n <k where k := 2+ ■ It is easy to see that A + intersects the top 
half- meridian according to the following ordered sequence of points: yo, Pi , f 2 _ , i •••• Similarly, 
intersects the top half- meridian at yo, ffi -f^) ^3") •••• Moreover, let y^ q and be the antipodal points 
of the equator which are the respective first intersections of A + and A~ with the equator. Note that 
they are reached at the same time T eq . Finally, consider the open subset of the top-hemisphere bounded 
below by the equator and obtained by removing the supports of A + and A~ up to time T eqi i.e. all the 
L%. That set is the disconnected union of the two "snake-shaped" simply connected regions S + and S~ 
(defined so that each S £ contains the center of rotation of Xg). Clearly S + and S~ are made of open 
segments of the meridian and open simply connected regions D £ k C NI £ defined as follows. For k = 0, Dq 
is delimited by Supp^p) and the segment [Poj-Pf] and, for k > 1, D% is delimited by L%_ 1 on the top, 
L e k on the bottoms and by the segments [P k ~I 1 , P)7 £ ] and [P| , P k+1 ] on the sides, see Fig. 6 B. 

In the sequel, if A, B are two subsets of points of 5 e , we say that A is above B (or equivalently B is below 
A) if A C D% and B C D%, with k < k', for some e , £ . 

Step 4 The switching curves (SC for short), associated to the set of extremals given in (40) and (41), are 
defined as follows: they can be divided in two families, {C k ) and (C k ). If e = ±, 1 < fc < A^o — 1 and 
s e [0, 7r], then 

Cf(s) = e x s v ^e x s es yo , C s k+1 (s) = e x ^C k £ (s). 

The boundary points of C| are C|(0) = an d CK 71 ") = Pr+i- By using Proposition 9 and since v(s) > ir, 
s £ [0, 7r], the support of Cf. is contained in the subset of Clos(NI £ ), delimited by the half-circle centered 
on 0P(a(2k + l)j and passing through the points P k , P k+1 , and the segment of the meridian [P k , P% +1 ], 
see Fig 6 B. In particular, a (SC) with boundary points in the top-hemisphere is entirely contained in the 
top-hemisphere and the intersection of its support with the top-meridian reduces to its boundary points 
(see Lemma 3). 

We next describe the shape of the first (SC) intersecting the equator. By symmetry, we may assume 
£ = +. We claim that its intersection with the equator reduces to the point -P(f) = (1,0, 0) T . Indeed, 
by the switching rules established in Proposition 9, the (SC) intersecting the equator is contained in 
{y e S 2 : 2/2 < 0, y 3 > 0} U {y e S 2 : y 2 > 0, y 3 < 0}. Taking into account the regularity of the (SC) and 
the values of its boundary points, the claim is proved, see Fig. 7 A. 
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4.4 Computation of Ns(ce) 

In the previous section, we provided detailed informations about extremal trajectories and switching curves 
but we did not show that every extremal of (40) and (41) is in fact time optimal. Anyway, a rigorous derivation 
of Ns(a) is possible with the available knowledge of time optimal trajectories combined with the subsequent 
lemmas. 

Lemma 7 Every time optimal trajectory 7 starting at y intersects the equator at most once. 

Proof of Lemma 7 we argue by contradiction. There would exist two distinct points of the equator qi,q; 
so that 7(t») = q%, j(tf) = q/ and j\(t t ,tt) i s entirely contained in the (closed) bottom hemisphere. Let j S ing 
be the integral curve of Fs (contained in the equator) connecting qi to q/. Consider now the region of the 
bottom hemisphere bounded by "fsing and 7|(t <) t,)- Taking into account, first, the relative positions of Xg, Xg , 
Fs and Gs along the equator and, second, the sign of fs in the bottom hemisphere, one can check that 
T (-fsing) < r(7|( ti)t )). The argument is similar to that of [32] (see also [12]) and is based on the use of Stokes 
theorem. Since time optimal trajectories starting at yo do not contain a singular arc, it follows that 7 cannot 
be time optimal. We reached a contradiction. H 

Lemma 8 Every time optimal trajectory 7, starting at yo and remaining in IH , is the projection of a time 
optimal trajectory of (£) starting at Id. 

Proof of Lemma 8 From the definition of the Hopf fibration, every trajectory 7 of (£)s, starting at yo, 
associated to an admissible control u and staying in IH , is the projection of the trajectory 7 of (£) starting 
at Id with the same control u. In particular, 7 and 7 have same time duration. It is clear that, if 7 is time 
optimal, then 7 is also time optimal. U 

Lemma 9 Recall that S e C PH. With the notations above, pick any point y in the region S E and let ~/ y be a 
time optimal trajectory connecting the north pole yo to y. If s g]0, 7t[ is the time duration of the first bang arc 
and T{y) the total time duration ofjy, then 7y|( Sj T(j/)] is entirely contained in S e . 

Proof of Lemma 9 By the switching rules of Proposition 9, along every time optimal trajectory contained 
in IH E , the control must switch from e to —e, when arriving at a switching curve Cf. In addition, the time 
optimal trajectory switches from being an arc of circle (integral curve of X%) to another arc of circle of bigger 
radius (integral curve of Xg e ). After rectification of the flow of Xg, (i.e. the one entering the (SC) C|), then, 
by taking into account Remark 9, one gets the situation depicted in Fig. 7 B. 

By contradiction, we assume that there exists a time optimal trajectory 7 with time duration T and first 
bang arc time duration s < T such that 7 connects yo to y G S £ and 7|( s ,t] exits from S s . Let t' be the smallest 
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Figure 8: Proof of Lemma 9. Here to fix the ideas we set e = 



— , £•' = +, k — I 



time (in [0, T}) so that j\(t',T] is entirely contained in S £ . Then j(t') belongs to Supp(A^ T j) U Supp(A^ T j) 
(see step 3 of Section 4.3 for the definition of T eq ). If j(t') is on the (top)-meridian, then it has to switch so 
that the interior bang time duration is constant, equal to 7r. Therefore 7|(f i T] will never re-enter S + . We thus 
deduce that j(t') is not on the meridian and, with no loss of generality, we will assume that j(t') belongs to 
the (1-dim.) interior of some LjJ", k > 1. Now we make the following two claims: 

Claim 1 With the notations above, there exist t" < t' < t'" such that 

7|(t',t'») C d( C S e and 7|(t",t') C d( +1 C S~ £ , /or some e' e {+, -}, 
i.e. 7 passes (backward in time) from S £ to S~ £ at time t' by going "down". 

Proof of Claim 1: it is clear that there exists a neighborhood U of t' so that "flu is an integral curve of Xg £ . 
Thanks to Remark 9 and to the argument above, intersects Int(L^) transversally (see figure 8) in such a 
way that 7, run backward in time, goes from D~£ to D~£ +1 . Claim 1 is proved. □ 

Now, by definition of t', j(t') G ^4 e |[ 0! T c ,]- Let j a b be the restriction of A £ between y and 7(i')- Consider 
7, the concatenation of "f a b and 7|(c j t]- The conclusion of Lemma 9 will follow if one can show that the time 
duration T" of 7 is less than T, the time duration of 7. This, in turn, amounts to show that T' , the time 
duration of 7 a & is less than t', the time duration of 7|[o,t']- This is the object of the next Claim. 

Claim 2 With the notations above, we have T" < t' . 

Proof of Claim 2: The trajectory 7, run backward in time from t', is an Xg £ -integral curve until it hits 
a (SC) in some D £ L e NI £ , for some integer L > k + 1, at a point C£ (s), s e]0, 7r[. One can easily conclude 
that the only possibility is L = k + 1. By Claim 1, a time optimal trajectory can pass (backward in time) from 
S £ to S~ £ only by going down, i.e. by passing from some D £ k to D £ k+1 . Therefore, by an elementary counting 
argument, one gets 

t! = s + t s + (k + l)v{s), 
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where t g is the time needed to go from j(t') to C~£ +1 (s). On the other hand, 

T' = (K+ 1)tt - 1, 

where t €]0, 7r[ is the time needed to A £ to go from j(t') to -P fc f fl . Since v(s) > n, then T' < t'. The proof of 
Lemma 9 is finished. B 



Remark 10 Coupled with the proof of Claim 2, a simple continuity argument implies that A + and A 
time optimal trajectories in the top hemisphere. 

Gathering all the information on time optimal trajectories, we are now able to compute Ns(ct). 

Proposition 13 For a e]0,7r/4[, we have 

N s (a) := 2 



are 



" 7T " 






' 7T " 


7T " 


-8tt- 




2 


-8a- 


4a. 



(42) 

Proof of Proposition 13 Let y E 5 + and 7, a time-optimal trajectory connecting 2/0 to y. The point y belongs 
to some , k < N and, by Lemma 9, 7 remains in 5 + . Since the function v takes values in [ir,ir + 7r/2], 
it is easy to sec, from Eqs. (40) and (41) and Remark 10 that 7, run backward in time, will go through the 
ordered sequence of regions D~£ , D^_ lt D~^_ 2 , etc, until hitting one of the two curves 7^ or j~ . Moreover, 
in each of the region Df , 7 will switch exactly once, thanks to Proposition 9. Therefore, the number of times, 
where an optimal trajectory 7 starting at yo switches, is exactly equal to the number of times 7 crosses the 
subset of the meridian contained in Nl. The same conclusion holds for points belonging to S~ , A + and A~ . 
By a systematic examination of all the possibles cases, we end up with (42). Note that Ng(a) is the number 
of switchings for a time optimal trajectory ending on the equator. H 



4.5 Geometric Remarks 

4.5.1 Relations with the Linear Pendulum 

In a fixed neighborhood of the north pole, the control system on the sphere (26) behaves, when a > is small 
enough, as a controlled linear pendulum. More precisely, let us consider the stercographic projection of the 
sphere from the south pole (0, 0, —1) on V, the tangent plane to the sphere at the north pole. If 2/1,2/2, 2/3 are the 
coordinates of the three dimensional Euclidean space where the sphere is embedded, a system of coordinates 
on V is (2/1,2/2), see Fig. 9. Let x$ and x^ be the projections of the equilibrium points of Xg and Xg in NI. 
An alternative way of parametrizing this problem (instead of fixing the radius of the sphere and varying the 
axes of rotations) consists of fixing the points Xq = (1,0) T , Xq = (— 1,0) T and varying the radius r of the 
sphere. The relation between a and r is tan(a) = 1/r. The range a e]0,7r/4[ becomes r e]l,oo[ and a — > 
corresponds to r — > 00. In V, fix a ball B(0,r o ) of radius r > centered in the origin, and consider the 
stereographic projection of the integral curves of Xq and Xq . For r — > 00, they become circles centered at the 
points Xq. Then, one easily sees that, in 5(0, r ), the limit system (and the associated synthesis) corresponds 
to a controlled linear pendulum (with the associated synthesis) of equation: 2/1 = —2/2, Vi = 2/1 ~ u , \ u \ < 1- 
Notice that linic^oo v(s) — ir, that is exactly the time duration of interior bang arcs for the linear pendulum. 

4.5.2 The time optimal problem on SU(2) 

The optimal control problem on NH is the projection (by a Hopf fibration) of an optimal control problem on 
50(3). Similarly, the corresponding problem on the whole sphere 5 2 is the projection (by an appropriate Hopf 
fibration) of an optimal control problem on SU(2). Indeed, SU(2) is the universal (double) covering of 50(3) 
and they have the same Lie algebra so(3). The existence of that double covering justifies, by a factor two, the 
difference between our bound and the bound (2), on the maximal number of switchings for the control problem 
on 50(3). Indeed, the index theory developed by Agrachev and Gamkrclidzc in [3, 5] provides a bound on 
the number of switchings by proving that a certain extremal is not optimal because it loses local optimality 
working at the Lie algebraic level. This is why the upper bound in (2) corresponds (essentially) to a control 
problem on SU(2), and thus, after projection, on a control problem on the whole sphere 5 2 , and not just on 
NH . The other factor two, of the difference between our bound and the bound given in (2), comes from the 
fact that in [3] the index of the second variation was estimated up to an additive factor 1 (see [3], p. 275). 
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Figure 9: Stereographic projection and synthesis of the linear pendulum 



5 Conclusion and Open Problems 

In the previous Section, we derived a set of properties of the optimal synthesis that were sufficient to compute 
the maximum number of switchings of a time optimal trajectory joining yo to any point of the north hemisphere. 
This enabled us to provide a precise estimate for N(a), a g]0,7t/4[. However, the following questions remain 
unsolved: 

Question 1 arc all the extremal trajectories (40) and (41) optimal in the north hemisphere? 
The answer to this question depends on the answer to the next question: 

Question 1' in the north hemisphere, are the switching curves C£(s), s G]0, 7t[, locally optimal? (The points 
s — 0, s — 7r are not included since we already know that the two abnormal extremal are optimal in 
M.) 

Roughly speaking we say that a switching curve is locally optimal if it never "reflects" the trajectories. More 
precisely, we have the following definition (clarified by Fig. 10). 

Definition 7 Consider a smooth switching curve C between two smooth vector field Y\ and Y% on a smooth 
two dimensional manifold. Let C(s) be a smooth parametrization of C . We say that C is locally optimal if, for 
every s € Dom{C), we have 

C{s) 7^ a\Yi{C{s)) + U2Y2(C(s)), for every a±,a2 s.t. aia2 > 0. (43) 

The points of a switching curve on which relation (43) is not satisfied arc usually called "conjugate points" . 

Remark 11 Notice that, if all the switching curves are locally optimal in the north hemisphere, it follows that 
the set of extremals (40) and (41) (restricted to NT) is an optimal synthesis for the problem (26) on K.P 2 . In 
this case, on RP 2 , the extremals (40) and (41) lose global optimality before losing local optimality. 

Question 2 If the answer to Question 1' is yes, what about the same question for the optimal control problem 
on S 2 1 More precisely, one would like to understand how the extremal trajectories (40) and (41) are 
going to lose optimality in a neighborhood of the south pole (i.e. if the loss of optimality is local or just 
global). 

Question 3 What is the shape of the optimal synthesis in a neighborhood of the south pole? 
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Figure 10: Definition 7 



In this section, we present the results of some numerical simulations which provide some hints regarding the 
above questions. More precisely: 

• There is strong numerical evidence for a positive answer to Question V . This means that the switching 
curves in the north hemisphere never reflect trajectories. In other words, situations like those considered 
in the proof of Lemma 9, (cf. Fig. 8) are not possible. 

• As concerns Question 2, we conjecture the following: 

CI The curves C|(s), s g]0,tt[ are locally optimal if and only if X+(C e k (0)) = cnXgCffl)) and 
Xs (Cf.(ir)) = ot2X s ~C k { / n)) with a.\,a.i > but not both vanishing. 

This condition is verified if and only if: k < f 2 ^] — 1, which simply follows from Remark 9. 

Set Na '■= [ j-] ■ Analyzing the evolution of the minimum time wave front in a neighborhood of the 
south-pole, it is reasonable to conjecture that: 

C2 For T < (Na — the synthesis built above is optimal. Every x e S 2 is reached in time T < 
(Na+1)tt. Every optimal trajectory has at most Na switchings and there exists an optimal trajectory 
having Na — 1 switchings. 

In the top of Fig. 11, the optimal synthesis is plotted. 

• Regarding Question 3, numerical simulations suggest that the shape of the optimal synthesis for time 
T > (Na — l)7r, depends on the remainder: 



r := tt — 2uNa = tt — 2a 



TT 

~2a~ 



Notice that r belongs to the interval [0, 2a[. More precisely, we conjecture the following: 

C3 For a s]0, 7r/4[, there exist two positive numbers et\ and 012 such that < ol\ < a < 012 < 2a and: 

CASE A: r e]a2,2a[. The switching curve starting at glues to an overlap curve that passes 

through the origin (see the bottom of Fig. 11, Case A). 
CASE B: r € [cui, 02]. An overlap curve starts exactly at P£ and passes through the origin. 
CASE C: r <g]0, a\[. The situation is more complicated and it is depicted in the bottom of Fig. 11, 

Case C. 

For r = 0, the situation is the same as in CASE A, but for the switching curve starting at P£ _ v 
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Figure 11: The Time Optimal Synthesis on the Sphere (top) and Optimal Synthesis in a neighborhood of the 
south pole, Case A and Case C (bottom) 
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